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We study the effects of a thin gaseous accretion disk on the inspiral of a stellar-mass black hole 
into a supermassive black hole. We construct a phenomenological angular momentum transport 
equation that reproduces known disk effects. Disk torques modify the gravitational wave phase 
evolution to detectable levels with LISA for reasonable disk parameters. The Fourier transform of 
disk-modified waveforms acquires a correction with a different frequency trend than post-Newtonian 
vacuum terms. Such inspirals could be used to detect accretion disks with LISA and to probe their 
physical parameters. 
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The inspiral of a stellar- mass compact object (SCO), 
such as a black hole or neutron star, into a supermassive 
black hole (SMBH) is among the most interesting gravi- 
tational wave (GW) sources for the Laser Interferometer 
Space Antenna (LISA) Extreme mass ratio inspirals 
(EMRIs) can be produced throiigh fragmentation of ac- 
cretion disks into massive stars 0, Q or through capture 
of compact remnants by hydrodynamic drag which 
are believed to be mass-segregated in galactic nuclei Q , 
as well as through other channels @ . Stars which reside 
within an accretion disk will lead to EMRIs, provided 
they become a SCO in less time than their inward mi- 
gration time. Although the expected EMRI event rate 
is rather uncertain (between a few tens to hundreds over 
EISA's lifetime, including coalescences and inspiral-only 
events 01)5 8- detectable fraction may originate in active 
galactic nuclei (AGNs) with an accretion disk. 

Accretion disks are efficient at extracting orbital an- 
gular momentum from the extreme mass-ratio binary. 
The SCO torques the disk gravitationally, inducing spiral 
density waves that carry away angular momentum Q. 
In planetary disks, the same phenomenon leads to mi- 
gration of planets towards their parent star. Planetary 
migration has been classified into different types (deter- 
mined by disk parameters, the EMRI's mass ratio and 
orbital separation) to distinguish circumstances where a 
gap opens around the planet (Type-II) from those with- 
out a gap (Type-I). In EMRIs, migration becomes the 
dominant source of angular momentum transport at sep- 
arations > 100A/„ where M, is the SMBH mass 0, Q 
(we use units G = c = 1). 

Migration changes the relation between the binary's 
binding energy and the GW luminosity, and hence it af- 
fects the inspiral rate and the GW phase evolution. EM- 
RIs enter the LISA sensitivity band only inside < 50M,, 
where GW angular momentum transport is dominant. 
Thus, migration acts perturbatively in LISA EMRIs. In 
this Letter, we examine whether the imprint of migration 
on the EMRI GW observables is detectable by LISA. In 



a companion paper [10[, we consider a broader range of 
disk effects and their impact on GWs in more detail. 

Disk Properties and Migration. — We consider 
radiatively-efficient, geometrically thin accretion disks, 
whose two most important free parameters are the accre- 
tion rate M, (overhead dots denote time-derivatives) and 
the a- viscosity parameter. AGN observations suggest 
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an accretion rate M, = m.M.Edd G (0.1, l)M,Edd 
Evidence for the magnitude of a is inconclusive, 
with plausible theoretical and observed ranges in 
(0.01, 1) [llj. We focus on Shakura-Sunyaev a-disks 



and /3-disks [1J|, which differ in whether viscosity is 



proportional to the total pressure (gas plus radiation) 
or only the gas pressure, respectively. This affects the 
surface density (E oc r^/^ and r~^^^ for a- and /3-disks 
when opacity is dominated by electron scattering). 
The local disk mass is much larger for /3~disks at radii 
r ^ 10'^ M,, leading to a stronger GW imprint. 

In the absence of a gap, Type-I migration models for 
angular momentum transport have been formulated [15, 
[Tij l but they are very sensitive to opacity and radiation 
processes [17| and lack the stochastic features observed 
in magnetohydrodynamic simulations jlq . The presence 
of a gap leads to Type-II models for angular momentum 
exchange T^l^]. These also oversimplify the process, as- 
suming either a steady state or quasi-stationarity. Type- 
II migration can also cease interior to a decoupling radius, 
Td , in the late stages of the inspiral, when the gas accre- 
tion velocity outside the gap becomes slower than the 
SCO's GW-driven inspiral velocity [lH. Alternatively, 
the gap can refill by non-axisymmetric or 3D inflow, 
restoring viscous torque balance from inside and outside 
the SCO's orbit and slowing the gaseous migration [2^ . 
Migration is mostly unexplored in the regime relevant 
to LISA EMRIs, i.e. for radiation-pressure dominated, 
optically-thick, geometrically-thin, relativistic, magne- 
tized and turbulent disks, with the SCO's mass rUi, ex- 
ceeding the local disk mass. 

Astrophysical uncertainties regarding accretion disks 
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TABLE I. Disk parameters for Type I and II migration models 
in 01 and /? disks. 



and migration in the regime relevant to LISA EMRIs 
lead us to consider a general power-law relation. 
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where is the SCO's rate of change of specific angu- 
lar momentum, ^gw is the loss due to GWs, and 5i a 
correction induced by migration. The power-law form in 
the reduced radius f = r/M, involves an amplitude A, 
which is parameterized in terms of normalized accretion 
disk (ai = a/0.1, rh, = m,/0.1), and mass parame- 
ters [M,5 = M./(105Mo), m^i = 771^/(10 Mq)]. The 
power-law indices (A^^Q: B) are given in Table U for rep- 
resentative migration models: rows I~2 correspond to 
Type I [1^ , 2-4 to steady state Type-II [l^ , 5 to quasis- 
tationary Type-II migration in the asymptotic limit for 
small f l20j (the latter is available for /3~disks only) . The 
gap decouples and Type-II migration ceases {A w 0) in- 
terior to fd = 1.4 X lQ~^a'^'^rn~^M~^m1^\^ for a and 
\ba^^^^^m^i^^^M~^^'^'^ml/^'^X^/^^ ford disks, where Xr 
is the gap radius (we adopt A = 1.7 [23|). Since disk ef- 
fects become stronger at larger radii, B > 0. 

GW Implications. — The change in the angular mo- 
mentum dissipation rate due to migration modifies the 
GW evolution, leading to a change in the accumu- 
lated GW phase and spectrum. For circular orbits, the 
quadrupolar GW phase can be computed from (/i>Gw — 
2 dr^t~^ d£/dr, where the orbital frequency is 17 ~ 

(Af./r'^)^/^, the binary's specific angular momentum is 
£ = r^Jl = My^r^/^, while the specific angular mo- 
mentum flux is given by Eq. ([T]). For a fixed final 
EMRI separation rf and observation time Tobs, the initial 
separation Tq is different from what it would be in vac- 
uum, as the radial inspiral evolution f is determined by 
4: r = J^/ L{d£/dr)-^dr. For an unperturbed EMRI, 
(f//fo)-4 « 1 + 33 (m,i/M25)(Tobs/yr)(f//10)-4. 

The correction to the GW phase, S(j)cw = 4>gw ~ 
where ^Jw is the accumulated phase in vacuum, is then 

(3) 
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where ^ = -(3 x 4-1/2 ^ 55)1^4 ^ B)-^{5 + 2B)-^A, 
X = rf/rQ, fo = tq/M, and we have expanded in S£ ^ 
1, which holds in the LISA regime. For fixed Tobs, we 
find that |(5(/)gw| increases and decreases with m^, for the 
Type-I and II models of Table |T1 respectively. 
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FIG. 1. Top: contours of logjQ |5(;/)gw| for different interaction 
models of two EMRI systems observed for one year. Bottom: 
liJt/iGwl for different f/ for the 5 models in Table H] Many of 
the models - especially those resembling /3-disks and Type II 
migration - shown in Table |I] (marked with symbols in the 
top, and individual curves on the bottom) lie well above the 
LISA sensitivity level (thick magenta lines). 

The top panel of Fig. [T] plots the dephasing in Eq. ([3]) 
for two typical LISA EMRIs at fixed r/ > rd as contours 
for different torque parameters (A, _B). The specific mi- 
gration models defined in TableUwith ai = 1 = m,i are 
marked with symbols. The bottom panel shows (J^gw 
for those models but with different ry, fixing Tobs = lyr 
(c.f. LISA'S planned lifetime is 3 years). For comparison, 
we also plot the total GW phase accumulation [O(IO^) 
top, thin line] and a rough measure of EISA's accuracy to 
phase measurements: (50gw > 10/p, where p is the sig- 
nal to noise ratio p{K) = 4 d,r{dj j dr)\h\^ S;,\] {r% 

with S'„[/(r)] the LISA detector noise [IJ] and h the 
Fourier transform of the orientation-averaged GW sig- 
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nal. We evaluate p at 1 Gpc (or redshift z 0.2; thick 
solid line) and at 10 Mpc (or z w 0.002; thick dashed 
line). For p < 10, we assume the EMRI is not detected 
at all, which explains the sharp rise in the detection level 
beyond a certain r/. Migration with a gap (empty sym- 
bols) causes a bigger phase shift because of the pileup 
of mass outside the gap. For < f/ < 50 but fixed 
{A, B, M,, m^, Tobs), the phase shift is constant within a 
factor ^ 3, but it quickly drops off for the Type-II models 
interior to the gap decoupling radius where A ^ 0. 

The Newtonian estimates presented here suggest 
that LISA EMRI observations might be able to probe 
accretion-disk induced migration. Figure [T] shows that a 
large sector of parameter space (^4, B) exists where the 
dcphasing is large enough to be detectable, and 6(j)Gw 
is very sensitive to the disk model and its parameters. 
One might worry, however, that the estimates in Fig. [T] 
are inaccurate due to the use of a Newtonian waveform 
model. We have verified that this is not the case through 
a relativistic waveform model that employs the calibrated 
effective-one-body scheme [25|. We have generated 1 
year-long waveforms for the systems plotted in Fig.[T]and 
included modifications to the radiation-reaction force due 
to migration, as parameterized by Eq. ([2]). Overall, we 
find the Newtonian results to be representative of the 
fully relativistic ones [loj . 

Just because migration produces a sufficiently large 
phase correction does not necessarily imply that LISA 
can measure it. For that to be possible, migration phase 
corrections must be non-degenerate, or at worst, weakly 
correlated with other system parameters. One can study 
if this is the case by computing the Fourier transform 
of the GW observable. We employ the stationary phase 
approximation (SPA) [2^, where one assumes the GW 
phase varies much more rapidly than the amplitude. The 
Fourier transform of h{t) = A{t) {t)] can then be 

approximated as h{f) = A{f)cxp[iil!{f)], where A{f) = 
(4/5)A[t(/)]/-i/2 and V(/) = 2^/<o - 0Gw(io), where 
/ is the GW frequency and to is the stationary point, 
defined by 27r/ = (d(/)Gw/rfi)t=to M- 

The corrections due to migration on the Fourier trans- 
form of the GW phase in the SPA, S^p = tp ~ i^vac, are 
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where we have defined A = -22-8^/551-85/5(^4 _^ 
i3)-^(5 + 2i3)-^Aexp(6.46i?), the normalized symmetric 
mass ratio f\ = m^i/Af.s and u = (7rA^/)/(6.15 x 10~^), 
and where A4 = ni'l/^Mt^^ is the chirp mass and 
V-vaT* = (3/128)(7rX/)-5/3 is the leading-order (New- 
tonian) vacuum Fourier phase. The amplitude of the 
SPA Fourier transform is corrected in a similar fashion: 

5|/l|/|/l|vac ~ <5V'/V' vac ■ 

Equation ^ is to be compared with the intrinsic gen- 
eral relativity (GR) corrections to the vacuum Fourier 



GW phase: ■0vac/'0- 
aq{mi,, M,), and the modulation induced by the orbital 
motion of LISA around the Sun. Migration corrections 
lead to negative frequency exponents in the Fourier phase 
(in Eq. (g]), -2S/3 < 0), while GR, post-Newtonian cor- 
rections in vacuum lead to positive powers of frequency, 
while the detector orbit is periodic with a 1 yr period. 
For a sufficiently strong signal, this suggests it might be 
possible to separate the migration effects from the other 
GR and detector orbit induced phase corrections. 

Modified gravity theories might introduce corrections 
similar to Eq. (|4]). The parameterized post-Einsteinian 
(ppE) framework [23|, devised to search for generic GR 
deviations in GW data, postulated such a phase modi- 
fication, allowing for both positive and negative B. De- 
generacies between disk and modified gravity effects with 
negative frequency exponents {B > 0) could then exist 
(e.g. Brans-Dicke theory or G{t) theories). The latter, 
however, have already been greatly constrained by bi- 
nary pulsar observations [l^. Moreover, alternative the- 
ory modifications should be present in all EMRIs, while 
disk effects will be present in only a small subset. 

A precise measure of whether a migration-modified 
waveform hi{f) is distinguishable from a vacuum wave- 
form /i2(/;A), where A stands for all disk parameters, 
requires a detailed Monte Carlo study that maps the 
likelihood surface. A rough measure of distinguishability 
can be obtained by calculating the signal-to-noise ratio 
(SNR) of the difference between a vacuum and a non- 
vacuum waveform p{5h) by minimizing only over a time 
and a phase shift. Using this crude measure, we demon- 
strate in the bottom panel of Fig. [5] that most of the 
migration models of Table Ulead to p{5h) > 10 within 4 
months of observation for a source at 1 Gpc. 



10" 

10' 
10° 
10' 

10-^ 
10-^ 

10' 
10° 
10' 
10-^ 
10-^ 
10^ 



^ SNR=2.9 SNR=5.9 SNR=8.6 SNR=11 



SNR=14 SNR=17 



M 1 ' 1 ' 1 ' 1 ' 1 
_! I l__ 

;r 1 1 




1 


I'll! 

■■■■■■■;V^i^^;|i 

n 
1 

1 ! 

\-. 


! 1 1 

■ '1/ ^ 








^--^^--^'-'f --JM-r^Sfl 


%/.■ 1 1 

m ^ 





1 = 

In 

1 : 

1 ! 

■ 



I (A,B)= 
( (A.B)= 
, (A,B)= 
, (A,B)= 
, (A,B)= 
. (A,B)= 
- (A.B)= 
(A.B)= 
(A.B)= 
(A,B)= 
(A,B)= 
(A.B)= 



le-12,5) 
le-13,5) 
le-1 1,4) 
le-12,4) 
le-13,6) 
le-1 4,6) 
Ib-14,7) 
le-16,7) 
le-16,8) 
le-17,8) 
le-17,9) 
le-18,9) 



I 



3 4 5 6 7 8 9 10 11 12^ 
t [Months] 



: (A.B)= 
f (A.B)= 
* (A,B)= 
» (A,B)= 
t (A.B)= 
3 (A,B)= 
3 <A,B)= 
i <A.B)= 
I <A,B)= 
I (A,B)= 



i-9,3) 
rlO,3) 
;-8,2) 
r9,2) 

!-7,l) 

^-10,4) 
^-6,25/8) 
■.-1,111) 
^-19,8) 
!- 13,5.9) 



FIG. 2. p{5h) as a function of observation time. Observe that 
p{Sh) > 5 within 1 year for a large set of {A, B). 

Going beyond TableU there exists a large sector of disk 
parameter space {A, B) for which the SNR of the wave- 
form difference exceeds threshold p{5h) > 10. Figure [2] 
plots p{Sh) at 1 Gpc as a function of observation time 
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for different values of {A,B) in Eq. ([2]). We also indicate 
(with labels over the vertical dashed lines) the SNR of the 
vacuum waveform at 1 Gpc and T = (0.5, 2, 4, 6, 9, 11.5) 
months. We calculate the waveforms with the relativis- 
tic model of ^ and (M.5,m*i) = (1,1), SMBH spin 
angular momentum \S,\ = 0.9M^ coaligned with the or- 
bital angular momentum, and initial and final separa- 
tions {ro,rf) ~ (24.5, 16)Af,. Observe that for a large set 
of disk parameters {A,B), p{Sh) > 10 within a one-year 
observation. Fitting to the smallest A with p{Sh) > 10 
for fixed B, we find that for these masses and orbital 
radii, LISA could measure logj^pA > ai + a2B, with 
ai = -5.7 ± 0.4 and aa = -1.4 ± 0.2. 

Discussion. — The GW observation of EMRI signals 
with LISA could be used to probe the uncertain physics 
of accretion disks. In particular, spiral density waves gen- 
erated by an orbiting SCO can transfer sufficient orbital 
angular momentum to alter the GW signal at levels that 
are detectable by LISA. The effects are strongest for pa- 
rameter choices resembling Type II migration (i.e. with 
the SCO opening a gap) in relatively massive (3 disks. A 
very crude (diagonal) Fisher analysis suggests that LISA 
could measure certain sectors of disk parameter space to 
better than 10%, for vacuum SNRs larger than 10, the 
details of which will be presented in a companion pa- 
per [lol | . This is no surprise considering that (5(/)gw is at 
worst ^ 10 times higher than LISA's sensitive curve in 
Fig. [TJ Detection of the predicted migration effect would 
reduce the uncertainty in existing theoretical models and 
offer the potential for extending the discussion to more 
complicated geometries (such as EMRIs with eccentric 
and/or inclined orbits). 

The detection of EMRIs in AGNs and the extraction 
of disk-parameters improve the prospects for finding elec- 
tromagnetic counterparts in the LISA error volume with 
consistent luminosities [l^. Coincident measurements 
would also allow EMRIs to serve as standard sirens to 
independently test cosmological models js^]. LISA EM- 
RIs are low-redshift events, for which weak lensing errors, 
dominant in comparable mass, SMBH standard sirens 
at highcr-z, are subdominant [29|. Disk effects will not 
compromise the ability to constrain cosmological param- 
eters, as they enter the GW observable with a differ- 
ent frequency signature, and are thus weakly correlated. 
Migration effects may also deplete the unresolved low- 
frequency EMRIs that contribute to the GW confusion 
noise background [23 |. 
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